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ABSTRACT 

We  present  a  new  lidar  method  for  investigation  of  the  aurorally-modified  ionosphere.  Resonance  nitrogen- 
ion  (N2+  XB)  lidar  system  is  capable  of  making  measurements  of  aurorally-ionized  nitrogen  in  the  polar 
ionosphere  up  to  altitudes  of  ~ 300  km.  We  present  a  description  of  the  prototype  resonance  lidar  system  that 
is  currently  being  developed  at  Poker  Flat  Research  Range.  We  compare  the  technical  features  and 
capabilities  of  this  lidar  system  to  conventional  resonance  lidar  systems  that  measure  atomic  metals.  Unlike 
incoherent  scatter  radar,  the  lidar  measures  a  specific  ionic  species  as  opposed  to  a  total  ion  profile.  We 
describe  how  assimilative  observational-modelling  studies  of  the  auroral  E-region  that  combine 
measurements  made  with  the  lidar  system,  a  radar  system  (i.e.,  the  Advanced  Modular  Incoherent  Scatter 
Radar),  and  a  meridian  scanning  photometer  with  a  multi-species  ionospheric  chemistry  model  will  be 
conducted.  We  discuss  how  these  observational-modelling  studies  will  yield  more  accurate  estimates  of  the 
ionic  populations.  We  discuss  how  these  studies  can  yield  estimates  of  nitric  oxide  production  while  it  is 
being  created  in  the  aurora,  estimates  of  the  auroral  particle  energy  spectrum,  and  allow  insights  into  the 
structure  of  pulsating  aurora. 


1.0  INTRODUCTION 

The  chemistry,  composition,  dynamics,  and  electrodynamics  of  the  ionosphere  (or  upper  mesosphere  and 
thermosphere)  are  significantly  influenced  by  the  aurora  [1,  2,  3].  Observations  of  the  aurorally  modified 
ionosphere  allow  insights  into  the  composition  and  circulation  of  this  region.  Auroral  electrons  and  protons 
are  accelerated  in  the  magnetosphere,  and  observations  of  their  low  altitude  impact  can  provide  significant 
insight  into  the  basic  neutral  and  ionic  composition  and  energetics  of  the  E-region.  Specifically,  auroral 
impact  on  thermospheric  nitrogen  creates  the  parent  N2+  ion  and  thus  feeds  subsequent  chemical  channels  to 
the  primary  E-region  ions,  NO+  and  02+.  Untangling  the  interplay  of  ionic  and  minor  neutral  composition  (i.e., 
NO)  in  this  chemical  cascade  remains  one  of  the  principle  challenges  in  auroral  aeronomy. 


Collins,  R.L.;  Su,  L.;  Lummerzheim,  D.;  Doe,  R.A.  (2006)  Investigating  the  Auroral  Thermosphere  with  N2+  Lidar.  In  Characterising 
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Nitric  oxide,  in  particular,  is  an  ill-determined  aurorally  generated  species  which  cannot  be  estimated  by  MSIS 
modelling.  Ground-based  estimates  for  its  abundance  and  altitude  distribution  can  more  fully  address  its 
lifetime  and  role  in  conversion  of  N2+  to  major  ions.  Synoptic  studies  of  NO  morphology,  based  on  dayside 
satellite  EUV  measurements,  have  concluded  that  increases  in  NO  at  mid-  and  high-latitudes  are  associated 
with  auroral  activity  [4]. 

Unusual  brightness  profiles  in  aurora  have  also  been  reported  but  remain  unexplained.  These  include 
extremely  thin  layers  (2  km)  of  emissions  in  pulsating  aurora  [5,  6]  and  enhanced  layers  near  the  lower 
borders  of  auroral  curtains  [7].  Determinations  of  the  altitude  and  altitude  extent  of  these  sharp  gradients  in 
the  auroral  brightness  have  heretofore  been  made  with  stereoscopic  optical  imaging. 

A  plethora  of  observational  methods  exist  to  make  auroral  observations,  ranging  from  satellites  or  in-situ 
measurements  by  sounding  rocket  to  ground  based  remote  sensing  using  radar  and  optical  methods.  Most 
observations  provide  line-of-sight  integrated  quantities,  which  must  be  supplemented  by  models  to  infer  the 
altitude  structure  of  the  aurorally  modified  thermosphere.  Tomographic  inversions  have  been  applied  to 
ground  based  and  satellite  based  optical  observations  to  reconstruct  the  three  dimensional  structure  with 
mixed  success  [8,  9,  10,  11,  12,  13].  Such  reconstructions  attempt  to  replicate  the  measurement  domain  of 
truly  range-resolved  diagnostics  such  as  incoherent  scatter  radar  (ISR)  and  lidars. 

ISRs  have  been  employed  for  over  30  years  to  measure  ion  and  electron  concentration,  plasma  temperature, 
and  ion  drift  profiles  in  the  thermosphere  [14].  Altitude  profiles  of  radar  measured  plasma  density  in  the 
auroral  E-region  can  be  related  with  few  assumptions  to  the  ionization  rate  profile  and  primary  electron 
energy  spectra  [15].  Although  the  brightness  ratio  of  auroral  emission  features  can  be  interpreted  in  terms  of 
mean  energy  of  the  auroral  electrons,  the  knowledge  of  the  altitude  profile  of  the  excitation  or  ionization  rate 
removes  the  uncertainty  of  an  assumed  energy  distribution.  Lanchester  and  co-workers  used  this  method  to 
conclude  that  the  very  intense  small  scale  structures  that  are  embedded  in  the  aurora  tend  to  be  better 
represented  by  narrow  Gaussian  energy  distributions  while  the  surrounding  less  structured  auroral  curtain  is 
best  represented  by  broader  Maxwellian  energy  distributions  [16]. 

While  ISRs  unambiguously  determine  the  electron  concentration,  the  technique  does  not  allow  unambiguous 
determination  of  the  species  of  ion  present  [6].  In  order  to  advance  our  understanding  for  the  cascade  of  N2+  to 
the  primary  E-region  ions,  and  to  characterize  the  interaction  of  ions  with  NO,  we  also  require  some  gauge  of 
ion  composition. 

Lidar  techniques  offer  the  ability  to  precisely  measure  distinct  molecular  species  and  to  differentiate  between 
the  different  vibrational  and  rotational  populations  that  are  present.  The  potential  of  lasers  to  probe 
thermospheric  molecular  species  has  been  recognized  since  the  1960s  [17].  Garner  and  Dao  [18]  presented  a 
more  detailed  analysis  of  the  problem  wherein  they  assessed  the  performance  of  a  lidar  system  based  on 
current  laser  technology  probing  a  1.5-km-thick  layer  of  aurorally-excited  molecules  and  ions  at  100  km. 
Collins  and  co-workers  [19]  extended  these  analyses  to  consider  profiling  the  auroral  species.  These 
researchers  surveyed  the  spectroscopy  of  molecular  nitrogen  and  determined  the  optimum  operational 
frequencies  for  the  lidar.  They  then  coupled  photometric  observations  of  the  aurora  with  an  ionospheric 
chemistry  model  and  determined  the  expected  profiles  of  several  aurorally-excited  species.  Collins  and  co¬ 
workers  used  this  analysis  to  establish  the  actual  requirements  for  an  optimum  auroral  lidar  system.  They 
found  that  the  optimum  lidar  system  was  a  resonance  nitrogen-ion  (N2+  XB)  lidar  system  that  would  excite 

nitrogen  ions  in  a  single-step,  N2+  (X2Z+g  (u  =  Ui))  — »  N2+  (B2Z+u  (u  =  u2)).  The  decay  and  remission  would 

occur  in  a  single  step  process  (the  "first  negative”  emission),  N2+  (B2Z+u  (u  =  o2))— »  N2+  (X2Z+g  (u  =  u3)).  A 
partial  term  diagram  of  molecular  nitrogen  is  shown  in  Figure  1.  (Note:  A  typo  in  [19]  incorrectly  labelled  the 

upper  state  of  N2+  as  B3fl+u).  Finally,  they  presented  a  simulated  set  of  observations  for  a  night  with  typical 
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auroral  activity  and  evaluated  the  performance  of  the 
resonance  lidar  system. 

Based  on  the  analysis  of  Collins  and  co-workers  [19], 
we  propose  to  conduct  an  assimilative  observational 
and  modelling  study  of  the  auroral  thermosphere  that 
combines  measurements  made  with  an  N2+  XB  lidar, 
an  ISR  (i.e.,  the  Advanced  Modular  Incoherent 
Scatter  Radar  (AMISR)),  and  a  meridian  scanning 
photometer  (MSP)  with  an  ionospheric  chemistry 
model.  We  show  will  show  in  Section  2  how  co¬ 
volume  measurement  of  total  ion  and  N2+  density  by 
ISR  and  lidar  techniques,  respectively,  constrain  a  set 
of  four  aeronomical  equations  to  predict  altitude 
profiles  for  major  ions  NO+  and  N2+,  as  well  as 
neutral  nitric  oxide.  In  Section  3.  we  discuss  the 
specific  implementation  of  lidar  and  ISR  techniques 
required  to  support  the  model  assimilation  and  predict 
the  performance  of  the  measurement  ensemble.  In 
Section  4  we  discuss  how  the  observations  and 
modelling  will  be  assimilated.  In  Section  5  we 
discuss  the  design  of  solid-state  lidar  systems  for 
studying  the  vibrational  population  of  nitrogen  ions  in 
the  auroral  thermosphere.  In  Section  6  we  present 
our  conclusions. 

2.0  FRAMEWORK  FOR  AN  ASSIMILATIVE  OBSERVATIONAL  AND 
MODELING  STUDY  OF  THE  AURORAL  THERMOSPHERE 

The  combination  of  lidar  and  radar  observations  provides  data  which  put  significant  constraints  on 
ionospheric  models.  To  demonstrate  the  interpretation  of  the  data,  assume  a  simplified  ionospheric  model  for 
the  auroral  E-region  based  on  the  most  important  reactions  that  are  listed  in  the  appendix  of  Rees  [2].  The  N2+ 
is  almost  entirely  produced  by  electron  impact  on  neutral  N2,  with  a  small  contribution  from  the  charge 

transfer  reaction:  +  N2  — »  02  +  N*  .  Neglecting  this  reaction  we  still  capture  90%  of  the  ionization  sources. 
The  dominant  loss  processes  are  reactions  of  N2+  with  molecular  and  atomic  oxygen.  The  importance  of  these 
reactions  depends  on  altitude.  A  very  small  fraction  (about  0.1%  in  the  auroral  E-region)  is  lost  by 
recombination.  This  leads  to  a  simple  continuity  equation  for  the  N2+  density.  The  ion-chemical  processes  are 
fast  compared  to  the  lidar  measurement,  and  assuming  steady  state,  the  continuity  equation  yields: 

'Tvt  =  [■ V2+]  (y5  [o2  ]  +  (/4  +  719  )E]) 

Where  PN2+  *s  the  auroral  ionization  rate,  square  brackets  denote  the  density  of  the  enclosed  species,  and  the 

reaction  coefficients  y  are  those  from  the  appendix  of  Rees  [2].  Measuring  the  N2+  density,  and  taking  the 
neutral  densities  from  the  MSIS  model  [20],  we  can  directly  derive  the  N2+  ionization  rate  as  a  function  of 
altitude.  The  brightness  of  the  N2+  first  negative  emission  at  427.8  nm  is  proportional  to  the  ionization  rate, 
and  the  integrated  emission  is  proportional  to  the  column  integrated  ionization  rate. 


Molecular  Nitrogen  Term  Diagram 


Figure  1:  Partial  term  diagram  for  molecular 
nitrogen.  The  transition  of  interest  is  the  first 
negative  system;  N2+  (X2Z+g  (v  =  vi))^N2+  (B2E+U  (v 
=  va))->  N2+  (X2I+g  (v  =  v3)). 
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The  E-region  ion-chemistry  for  the  major  ions  02+  and  NO+  is  more  complicated  but  can  also  be  simplified. 
02+  has  three  major  sources;  charge  transfer  from  N2+,  ionization  of  02  by  electron  impact,  and  charge  transfer 
from  0+.  The  latter  is  a  small  source  in  the  E-region  and  can  be  neglected.  The  mechanisms  by  which  02+  is 
lost  are  dissociative  recombination  and  reactions  with  NO  and  N2.  Assuming  steady  state  and  following  the 
same  notation  as  above  (ne  is  the  electron  density,  a  is  the  recombination  coefficient  from  [2])  we  get  the 
following  equation  for  the  02+  ionization  rate  from  the  continuity  equation: 

%.+rs[»2*][a]= 

(rl5[A'o]+  r,  ['V2]+al«,)[02  ] 

NO+  is  produced  only  by  chemical  reactions  involving  N2+,  02+,  0+,  the  major  neutrals  and  NO,  and  is  lost  by 
dissociative  recombination.  The  steady  state  continuity  equation  yields: 

T'*  ]  [C]+  y8  [OJ  ][iV2  ]  + 

r,s[0*][Af0]= 

a^ne^NO+\ 

In  addition  to  these  equations  we  have  quasi-neutrality: 

».-[jwr]+[o;]+M 

The  lidar  measures  the  N2+  density  from  which  we  derive  the  N2+  ionization  rate.  The  radar  measures  the  total 
electron  density  ne.  If  we  take  the  major  neutral  densities  from  MSIS,  we  then  have  three  equations  for  the 
four  remaining  unknowns,  [02+],  [NO+],  [NO],  and  P02+-  The  ionization  cross  sections  for  electron  impact  on 

02  and  N2  are  similar,  and  the  altitude  profiles  of  the  densities  of  these  two  neutral  species  are  similar,  so  we 
may  assume  that  the  02  ionization  rate  is  proportional  to  the  N2  ionization  rate.  The  proportionality  factor  can 
be  obtained  from  our  auroral  electron  transport  model  and  is  nearly  independent  of  altitude,  which  gives  us  the 
required  fourth  equation. 

Thus  combining  measurements  from  the  lidar  (N2+),  radar  (ne),  and  photometer  (column  integrated  N2+ 
emission)  with  several  assumptions  (steady  state,  and  negligible  0+  density  in  the  E-region),  we  can  solve  this 
simplified  E-region  model  and  obtain  the  altitude  profiles  of  the  densities  of  the  major  ions  (N2+  02+  NO+)  and 
NO. 

3.0  TECHNICAL  APPROACH 
3.1  Proposed  N2+  XB  Lidar  System 

Since  1999  an  excimer-pumped  dye  laser  has  been  employed  as  a  resonance  lidar  transmitter  at  Poker  Flat 
Research  Range  (PFRR),  Chatanika,  Alaska  (65°N,  147°W).  This  laser  system  has  been  employed  to  conduct 
broadband  measurements  of  the  mesospheric  sodium  (Na)  and  iron  (Fe)  layers  with  1  pm  spectral  resolution 
provided  by  the  manufacturer  [21].  The  dye  laser  oscillator  is  a  classic  Littrow  cavity  that  employs  a  rotating 
diffraction  grating  as  the  tuning  element.  In  2002  the  laser  grating  control  software  system  was  upgraded  to 
allow  higher  spectral  resolution  (0.1  pm)  measurements  of  the  Fe  layer  under  operator  control  and  make  initial 
Fe  Boltzmann  lidar  temperature  measurements  [22].  In  2004  a  microprocessor-based  closed-loop  tuning 
control  system  was  developed  to  allow  automated  high-resolution  tuning  of  the  laser  [23].  Progressive 
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improvement  of  the  lidar  system  has  allowed  researchers  to  conduct  multi-year  studies  of  the  mesospheric 
metal  layers  with  increasing  accuracy  and  precision  [24].  These  improvements  in  the  tuning  control  of  the  dye 
laser  are  critical  for  studying  species  that  are  variable  in  space  and  time  and  have  complex  molecular  spectra. 

To  date,  lidar  systems  at  PFRR  have  employed  telescopes  with  diameters  of  38  cm  and  60  cm.  In  2003  a  1-m 
class  (41  inches  single-mirror)  Cassegrain  telescope  was  transferred  from  the  University  of  Illinois  at  Urbana 
Champaign  to  the  University  of  Alaska  Fairbanks.  The  1-m  telescope  will  allow  measurements  with  signals  3- 
7  times  larger  than  previously  attainable  at  PFRR.  The  excimer-pumped  dye  laser  and  the  1-m  telescope  will 
serve  as  the  transmitter  and  receiver  respectively  of  this  new  N2+  XB  lidar  system.  Both  the  laser  transmitter 
and  telescope  receiver  are  proven  instruments  and  have  supported  a  variety  of  lidar  studies.  The  characteristics 
of  the  proposed  system  are  tabulated  in  Table  1.  The  transmission  and  reception  wavelengths  have  been 
chosen  based  on  the  spectroscopy  at  a  temperature  of  500  K.  As  the  rotational  distribution  of  the  molecules 
changes  with  temperature  the  optimum  lidar  wavelengths  change  by  several  nanometers.  The  ability  of  the 
excimer-pumped  dye  laser  to  tune  over  relatively  wide  wavelength  ranges  is  critical  for  this  prototype  lidar 
system  to  ensure  that  the  measurements  are  optimized  for  the  actual  geophysical  conditions. 

In  studies  of  the  mesospheric  Na  and  Fe  layers,  the  observations  are  facilitated  by  both  the  relative  stability  of 
these  mesospheric  species  and  the  well-defined  atomic  spectroscopy  of  Fe  and  Na.  Thus  lidar  measurements 
of  mesospheric  Na  and  Fe  can  be  optimized  using  the  resonance  echoes  from  the  mesosphere.  However,  the 
aurorally  produced  N2+  layer  is  more  variable  in  space  and  time  and  such  operational  methods  cannot  be  relied 
upon.  Furthermore,  the  standard  molecular  references  (i.e.,  Iodine  and  Tellurium  [25,  26])  do  not  cover  the 
390  nm  region.  The  N2+  XB  lidar  system  will  use  a  standard  references  (e.g.,  stabilized  HeNe  laser)  and  a 
tunable  Fabry  Perot  Interferometer  as  a  tuning  reference  for  the  laser  transmitter. 

In  studies  of  the  mesospheric  Na  and  Fe  layers,  where  the  transmission  and  reception  wavelengths  are  the 
same  (i.e.,  589  nm  and  372  nm  respectively)  the  concentration  of  the  metals  is  determined  by  taking  the  ratio 
of  the  resonance  signal  from  the  mesospheric  metal  layers  to  the  Rayleigh  signal  from  the  stratosphere.  This 
ratio  makes  the  resonance  lidar  measurements  more  robust  by  allowing  normalization  of  the  laser  power  and 
atmospheric  transmission  in  the  lidar  equation  [27].  The  transmission  and  reception  resonance  wavelengths  in 
the  N2+  XB  lidar  system  are  separated  by  -  lnm.  Thus  the  wavelength  of  the  Rayleigh  and  resonance  scatter 
is  separated  by  ~  lnm.  The  N2+  XB  lidar  system  will  employ  either  a  broadband  receiver  (>  lnm),  a  beam 
splitter  in  the  receiver  telescope  with  narrowband  filters  (~  0.3  nm)  for  simultaneous  detection  of  each 
wavelength,  or  a  cycling  approach  where  the  lidar  signal  is  recorded  at  one  wavelength  for  a  given  time 
interval  and  then  recorded  at  the  other  wavelength  for  a  given  time  interval. 


TABLE  1.  Proposed  N2+  XB  Lidar  System 


Transmitter 

Receiver 

Laser 

Excimer-pumped  dye 

Telescope  Diameter 

1  m 

Pulse  Energy 

50  mJ 

Detector 

Photomultiplier  Tube 

Pulse  Repetition  Rate 

20  pps 

Filter  Bandwidth 

0.3  nm 

Wavelength 

390.303  nm 

Wavelength 

391.537  nm 

Measurement  Resolution 

Temporal  Resolution 

300  s 

Spatial  Resolution 

5  km 
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3.2  Expected  Lidar  Performance 

Given  the  geophysical  variability  of  the  aurora,  the  performance  of  a  given  lidar  system  is  highly  dependent 
on  the  level  of  auroral  activity.  A  more  active  aurora  produces  higher  concentrations  of  excited  species  at 
lower  altitudes  and  yields  higher-quality  lidar  measurements,  while  less  active  aurora  yield  lower-quality  lidar 
measurements.  Collins  and  co-workers  [19]  reviewed  two  years  of  MSP  observations  and  selected  the  night 
of  6-7  February  1994  (UT  day  38)  as  representative  of  moderate  auroral  activity  with  varying  particle 
characteristics.  Figure  2  shows  the  brightness  of  two  emission  features,  the  “6300  line”  and  the  “4278  band”, 
over  the  course  of  the  night.  The  intensity  is  measured  in  Rayleighs  (R),  which  is  used  as  a  measure  of 
brightness  of  a  uniform  source,  such  as  the  airglow  or  aurora  [e.g.,  3].  The  aurora  on  the  selected  night  is 
brightest  early  in  the  night  (1900-2100)  and  just  after  midnight  (0000-0200)  when  several  discrete  auroral  arcs 
pass  overhead  from  north  to  south  and  the  intensity  in  the  green  line  (0(1S)^  0(1D)  557.7  nm)  reaches  about 
100  kR.  By  comparison,  early  in  the  night  (2100-2200)  the  green  line  intensity  is  about  1  kR,  which  is  at  the 
threshold  of  visual  perception.  One  MSP  channel  (486.1  nm)  monitors  the  Hp  brightness,  which  indicates 
proton  precipitation.  The  observed  Hp  brightness  on  the  night  of  6-7  February  1994  was  insignificant, 
indicating  that  the  auroral  brightness  was  exclusively  due  to  electron  precipitation.  The  average  N2  (A3E+  u) 
and  N2+  (X2S+g)  state  concentration  and  temperature  profiles  are  plotted  in  Figure  3.  The  peak  concentrations 
and  column  abundances  are  largest  (>109m'3  and  ~1015m'2  respectively)  when  discrete  auroral  arcs  pass 
overhead  (1930,  0039-0200  LST).  These  large  abundances  correspond  to  those  times  when  the  precipitating 
electrons  penetrate  lowest  into  the  atmosphere,  and  the  peak  of  the  ion  layers  are  at  their  lowest  altitudes. 
Recent  satellite  observations  suggest  that  the  topside  N2+  concentrations  may  be  significantly  higher  in  than 
these  model  results  indicate  [28]. 

The  expected  lidar  signal  was  calculated  every  300s  from  1902  LST  to  0657  LST  yielding  144  lidar  signal 
profiles.  The  lidar  signal  profile  with  the  largest  peak  amplitude  at  0057  LST  is  plotted  with  the  corresponding 
signal-to-noise-ratio  (SNR)  in  Figure  4.  The  signal  with  the  highest  SNR  occurs  at  0057  LST  during  the 
passage  of  discrete  auroral  arcs  overhead.  The  average  signal  and  SNR  profiles  for  the  whole  night  are  plotted 
in  Figure  5.  The  expected  lidar  signal  profile  has  an  SNR  greater  than  10  over  a  range  from  approximately  110 
to  200  km  for  these  300s-5  km  measurements. 


Figure  2:  Meridian  scanning  photometer  measured 
intensities  along  the  magnetic  zenith  at  630.0  nm  and 
427.8  nm  at  Poker  Flat  on  6-7  February  1994.  The 
measurements  are  made  at  16  s  intervals  with  a  field- 
of-view  of  1°. 


Temperature  (K) 

100  1000 


Figure  3:  Average  concentration  profiles  of  N2(A3£+U), 
N2+  (X2S+g),  and  temperature  for  6-7  February  1994. 
The  temperature  profile  is  taken  from  the  MSIS  model 
and  used  in  deriving  the  concentration  profiles  from 
the  MSP  data  in  Figure  2. 
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Figure  4:  Expected  photon  count  profile  (solid)  and 
signal-to-noise  ratio  (SNR;  dashed)  measured  by 
proposed  auroral  lidar  system  during  passage  of 
auroral  arcs  at  0057  LST  on  7  February  1994.  The 
photon  count  profile  represents  300s  of  integrated 
lidar  echoes  measured  at  5  km  resolution. 


SNR 


20  15  10  5  0 


Photon  Count 


Figure  5:  Expected  photon  count  profile  (solid)  and 
signal-to-noise  ratio  (SNR;  dashed)  measured  by 
proposed  auroral  lidar  system  averaged  over  whole 
night  of  6-7  February  1994.  The  photon  count  profile 
represents  300s  of  integrated  lidar  echoes  measured 
at  5  km. 


The  average  profiles  show  that  measurements  with  SNR  greater  than  10  occurred  between  110  and  190  km. 
The  altitude  distributions  of  the  measurement  occurrence  at  several  SNR  are  plotted  in  Figure  6. 
Measurements  with  an  SNR  greater  than  5  are  expected  to  be  made  >50%  of  the  time  at  altitudes  between  110 
and  260  km.  Measurements  with  an  SNR  greater  than  10  are  expected  >  50%  of  the  time  over  an  altitude 
range  of  120  and  180  km.  Figures  4-6  highlight  the  potential  of  resonance  lidar  for  making  species  specific 
measurements  in  the  E-  and  F-  Regions. 

The  analysis  is  based  on  lidar  data  acquired  at  a  constant  resolution  of  300s  and  5  km.  In  practice,  the  raw 
lidar  data  will  be  collected  at  a  higher  resolution  than  300s  and  5  km  to  allow  post-measurement  integration  at 
different  height  and  time  resolutions  to  yield  measurements  with  higher  lower  resolution  and  higher  SNR  or 
higher  SNR  and  lower  resolution,  in  order  to  investigate  different  and  distinct  features  in  the  thermosphere. 


3.3  Incoherent  Scatter  Radar  Measurements 

This  lidar  investigation  will  exploit  the  enhanced  beam 
pointing  capabilities  of  NSF's  planned  AMISR  facility  at 
PFRR.  The  AMISR  facility  establishes  a  new  state-of-the- 
art  for  IS  radar  design  by  implementing  fully  electronic 
beam  steering  with  a  phased  array  of  4096  UHF 
transceivers.  This  beam  pointing  capability,  available  on  a 
pulse-by-pulse  basis,  enables  the  formation  of  ionospheric 
plasma  'images'  with  temporal  resolution  limited  by 
available  signal  SNR  and  desired  precision  in  recovered 
plasma  density  and  temperature  estimates. 


In  order  to  address  the  science  goals  previously  outlined, 
the  AMISR  will  provide  contextual  measurements  for  the 
evolution  of  total  ion  concentration  and  plasma 
temperature  for  flux  tubes  threading  the  region  of  lidar 


Figure  6:  Expected  occurrence  frequency  of 
lidar  measurements  of  a  given  SNR  on  the  night 
of  6-7  February  1994.  Measurements  with  an 
SNR  greater  than  5  occur  over  50%  of  the  time 
over  the  altitude  range  of  120-270  km. 


RTO-MP-IST-056 


2-7 


UNCLASSIFIED/UNLIMITED 


Altitude  (km) 


UNCLASSIFIED/UNLIMITED 

Investigating  the  Auroral  Thermosphere  with  N2*  Lidar 


ORGANIZATION 


induced  N2+  emission.  Because  the  lidar  beam  is  pointed  into  the 
zenith,  as  shown  in  Figure  6,  N2+  backscatter  from  a  range  of 
altitudes  maps  to  a  spread  of  individual  field  lines.  The  AMISR 
platform  is  tipped  15  degrees  to  the  magnetic  north  and  is  designed 
to  provide  beam  steering  over  angles  ±30°  from  boresight  as 
indicated  with  gray  shading  in  Figure  7.  Thus,  plasma  parameters 
can  be  imaged  along  field  lines  that  intersect  the  lidar  beam  over  the 
100  to  300  km  range  of  altitudes  predicted  for  N2+  X2Z+g  species 
[19].  We  now  address  anticipated  AMISR  temporal  and  spatial 
sensitivity  within  this  cone  of  coverage. 

The  AMISR  has  an  effective  aperture  of  900  m2,  radiated  power  of 
2.0  MW,  a  system  temperature  of  120  K,  and  a  transmit  duty  cycle 
of  10%.  AMISR  performance  can  be  scaled  from  measurements  at 
the  Sondrestrom  IS  radar  facility  by  assuming  that  a  common  320 
ps  coded  pulse  mode  is  employed  to  provide  3 -km  resolution  in  the 
E-region  [29,  30].  Enhanced  AMISR  sensitivity,  relative  to 
Sondrestrom,  due  mainly  to  significantly  larger  antenna  aperture  (by 
a  factor  of  2.25)  and  transmitted  pulse  duty  cycle  (by  a  factor  of 
3.3),  is  moderated  somewhat  by  slightly  lower  radiated  power  (by  a 
factor  of  0.8)  and  a  slightly  higher  system  temperature  (120  K 
versus  85  K). 

The  relationship  between  desired  precision  in  plasma  parameters, 
SNR,  and  number  of  transmitted  pulses  has  been  prescribed  [31]. 
This  formula  was  applied  to  Sondrestrom  for  the  aforementioned  E- 
region  pulse  mode  with  assumption  that  plasma  parameters  should 
have  no  greater  than  10%  uncertainty.  This  relationship  was  then 
scaled  to  the  AMISR  system  attributes  and  the  transmitted  pulse 
number  converted  to  required  integration  time  with  knowledge  of 

the  interpulse  period  (29  ms).  Figure  8  provides  a  nomograph  for 
required  integration  time,  as  a  function  of  plasma  density  and 
altitude,  required  to  recover  AMISR  plasma  measurements  with 
10%  precision. 


3.4  Meridian  Scanning  Photometer  Measurements 

The  MSP  at  PFRR  simultaneously  records  the  airglow  in  four 
channels  (i.e.,  557.7  nm,  427.8  nm,  620.0  nm  and  486.1  nm)  [32]. 
The  MSP  scans  along  the  magnetic  meridian  from  horizon  to 
horizon,  completing  a  measurement  every  16  seconds  with  a 
instantaneous  field-of-view  of  1°.  Only  the  observations  in  the 
magnetic  zenith  are  used  as  these  represent  the  column  integrated 
emission  rates.  The  MSP  operates  in  an  ongoing  data-logging 
fashion  as  part  of  the  facility  instrument  suite  at  the  Davis  Science 
Center  at  PFRR. 
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Figure  7:  Viewing  geometry  of 
coordinated  N2+  XB  lidar  and  AMISR 
measurements  of  the  aurorally  modified 
thermosphere. 
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Figure  8:  Nomograph  for  required 
integration  time,  as  a  function  of  plasma 
density  and  altitude,  required  to  recover 
AMISR  plasma  measurements  with  10% 
precision. 
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4.0  COORDINATED  MEASUREMENTS  AND  MODEL  ASSIMILATION 

The  N2+  resonance  lidar  measurements  will  be  nested 
both  spatially  and  temporally  in  the  AMISR  and  MSP 
measurements.  The  MSP  provides  the  widest  view  of 
the  thermosphere  from  horizon-to-horizon  along  the 
magnetic  meridian  every  16s.  As  discussed  above,  the 
resolution  and  accuracy  of  the  lidar  and  radar 
Elevation  scans  of  total  plasma  density  measured  by 
the  ISR  at  the  Sondrestrom  Upper  Atmospheric 
Research  Facility  in  Kangerlussuaq,  Greenland  in 
January  2004  are  plotted  in  Figure  9.  The  four  panels 
show  successive  3 -minute  measurements  (top  left,  top 
right,  bottom  left,  bottom  right)  of  the  total  plasma 
density.  The  lower  two  panels  show  the  presence  of 
relatively  stable  ionization  with  a  density  of  4xl05 
cm"3  at  1 15  km.  From  the  nomograph  in  Figure  8,  we 
see  that  it  would  take  the  AMISR  system  2  seconds  to 
yield  a  measurement  profile  of  the  total  plasma 
density  with  10%  accuracy.  Assuming  that  the 
concentration  of  N2+  is  1%  of  the  total  plasma  density 
and  using  the  analysis  in  Section  3.2.,  we  expect  the 
lidar  to  yield  measurements  of  the  N2+  profile  with  a  resolution  of  1  km  and  30  seconds  with  10%  accuracy. 
Thus  for  the  radar  observations  in  Figure  9  we  envision  the  following  measurement  strategy;  during  a  single 
zenith-pointing  lidar  measurement  we  would  anticipate  15  AMISR  measurements  of  the  total  plasma  density 
profile.  The  AMISR  measurements  could  be  made  along  the  same  direction  as  the  lidar,  however  we  would 
propose  to  make  a  set  of  15  measurements  that  would  employ  the  rapid  electronic  steering  of  the  AMISR  to 
yield  the  evolution  of  the  spatial  structure  in  the  neighborhood  (i.e.,  ±30°  which  is  approximately  ±60  km  at 
115  km))  of  the  lidar  measurement.  During  the  same  30  second  period  the  MSP  would  yield  two  horizon-to- 
horizon  scans  of  the  airglow  along  the  magnetic  meridian  providing  the  spatial  and  temporal  evolution  of  the 
aurora  across  the  sky.  While  the  observation  methodology  of  the  MSP  is  fixed,  we  can  adapt  the  N2+ 
resonance  lidar  and  AMISR  measurements  to  the  specific  auroral  conditions.  Measurements  could  be  made  at 
higher  resolution  during  auroral  active  periods  to  address  studies  of  auroral  precipitation. 

The  coordinated  lidar,  radar  and  photometer  measurements  would  then  be  used  to  constrain  the  closed  set  of 
E-region  continuity  equations  (see  Section  2)  and  allow  determination  of  the  ionization  rates  and  altitude 
profiles  of  the  densities  of  the  major  ions  (N2+  02+  NO+)  and  NO.  The  development  of  a  computationally 
robust  and  physically  significant  assimilation  and  inversion  scheme  is  not  trivial.  This  development  will  draw 
on  analyses  of  the  model  simulations  and  a  suite  of  experimental  case  studies  that  represent  a  range  of 
geophysical  conditions  and  measurement  quality.  The  measurements  of  the  N2+  profile  can  also  be  used  in  ISR 
standard  inversion  methods  to  yield  estimates  of  the  primary  electron  energy  spectrum  [15]. 

As  we  have  noted  in  Section  3.2,  during  aurorally-active  periods  with  strong  auroral  precipitation  the  N2+ 
concentrations  increase  and  the  lidar  measurements  will  yield  N2+  concentration  profiles  at  higher  resolution 
(compare  SNR  in  Figures  4  and  5)  and  allow  studies  of  the  auroral  precipitation. 
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Figure  9:  Elevation  scans  of  total  plasma  density 
measured  by  the  ISR  at  Sondrestrom  on  18 
January  2004.  The  panels  represent  successive 
scans  taken  every  3  minutes. 
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5.0  MEASUREMENTS  OF  THE  DISTRIBUTION  OF  ROTATIONAL  STATES  IN 
THE  AURORA 

The  current  study,  employing  a  broadly  tunable  laser,  is  a  first  step  in  developing  resonance  lidar  systems  for 
investigating  the  thermosphere.  To  study  the  rotational  distribution  of  ions  in  the  thermosphere,  a  dual-laser 
lidar  system  could  make  simultaneous  measurements  of  two  rotational  states  and  study  the  temperature 
distribution  and  thermodynamic  equilibrium  of  the  auroral  thermosphere.  A  similar  system  has  been  employed 
to  measure  Boltzmann  temperatures  in  the  mesospheric  iron  layer  [33].  Such  a  system  would  employ  two 
tunable  vibronic  solid-state  lasers  (e.g.,  Alexandrite)  that  are  tuned  to  the  appropriate  lines  by  injection 
seeding  and  frequency  doubling  [34].  A  diagram  of  such  a  lidar  system  is  shown  in  Figure  10.  For  example 
one  laser  could  transmit  at  391.081  nm  (X2S+g  (J  =  0))  — »  B2S+D  (J  =  i))  while  the  other  laser  could  transmit  at 
390.303  nm  (X2S+g  (J  =  9))  — »  B2E+D  (J  =  i0)).  For  this  transmitter  configuration  the  detectors  would  be  tuned 
to  391.258  nm  (B2S+D  (J  =  1)  — >  X2S+g(J  =  2))  and  391.537  nm  (B2X+D  (J  -  iO)  ->  X2Z+g(J  =  11))  respectively. 
The  precise  choice  of  wavelengths  could  be  optimised  based  on  the  expected  rotational  distribution  of  the 
population  of  molecular  ions.  The  two  lasers  could  be  triggered  alternately  and  so  use  a  single  receiver.  Such 
an  approach  could  be  adapted  to  simplify  the  receiver  requirements  and  so  avoid  the  optical  losses  associated 
with  beam  splitters  and  the  cost  of  multiple  detectors.  Such  pulse-to-pulse  techniques  are  currently  employed 
in  sodium  Doppler  wind  temperature  lidar  systems  [35]. 


Rotational  State  Resonance  Lidar 


LIDAR  TRANSMITTER  LIDAR  RECEIVER 


Such  a  system  has  several  technical  advantages  over  the  prototype  system  with  a  single  transmitter.  Standard 
molecular  spectral  references  can  be  used  to  stabilize  the  continuous  wave  (CW)  seed  lasers  at  782  and  781 
nm.  Commercial  optical  materials  and  components  are  more  readily  available  at  780  nm  than  at  390  nm. 
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6.0  CONCLUSION 

In  summary  this  paper  presented  the  rational  for  the  development  of  molecular  resonance  lidar  systems  as 
tools  to  study  the  auroral  thermosphere.  These  lidars  provide  spectroscopic  measurements  of  specific 
molecular  ionic  species.  The  combination  of  lidar  and  incoherent  scatter  radar  measurements  and  assimilative 
models  has  the  potential  to  estimate  the  ion  composition,  nitric  oxide  and  ionisation  rates  in  the  thermosphere. 
A  prototype  N2+  XB  lidar  system  that  is  under  development  at  Poker  Flat  Research  Range  has  been  described 
and  compared  with  traditional  metal  resonance  (e.g.,  sodium)  lidar  systems.  This  prototype  lidar  system  will 
yield  measurements  of  at  a  resolution  of  300s  and  5  km  with  SNR  >  10  in  the  100-200  km  altitude  range  under 
moderate  auroral  activity.  More  active  aurora  will  yield  higher  resolution  measurements  over  a  larger  altitude 
range.  The  prototype  N2+  XB  lidar  system  will  provide  scientific  and  engineering  data  for  the  development  of 
dual-laser  lidar  systems  that  can  measure  the  distribution  of  rotational  states  in  aurorally  ionised  molecular 
nitrogen.  These  vibrational  state  resonance  lidar  will  allow  characterization  of  auroral  energy  deposition  in 
the  ionosphere. 
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Lidar  or  Laser  Radar 
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Receiver 


A  lidar  employs  a  laser 
transmitter  a  telescope  receiver 
and  associated  control  and 
acquisition  systems. 


UAF  students 
Tao  Hou  and 
Tim  Stern. 


Telescope 


Optics 


Detector 


Transmitter 


Reference 
&  Control 


Laser 


Computer 


Data 

^Acquisition  y 


Atmospheric  Lidar  Echoes 


Light  propagates  out 
from  the  laser  and  is 
scattered  back  into 
the  telescope. 

The  scattering 
mechanisms  include; 

Mie  Scatter 

Rayleigh  Scatter 

Resonance  Scatter 

Raman  Scatter 


14  November  1997 


102  103  104  105  106 


Photon  Count 

Rayleigh  and  resonance  lidar  signals 
plotted  as  a  function  of  altitude. 
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Lidar  Equation 

Number  of  photons  transmitted 


N(z)  =  (tiT) 


e<y 

A  )  +  NBRLAt(2Az/c) 

/ 

Background  Light  -  Noise 


4nz 


Efficiency 


Probability  of  photon  being  scattered 


r| 

T 

El 

Rl 

At 

h 


lidar  receiver  efficiency 
two  way  atmospheric  transmission 
energy  per  pulse  (J) 
laser  repetition  rate  (s'1) 
observation  interval  (s) 

Planck's  constant  (6.63  x  10'34  J/s) 


c 


’eff 


p(z) 

Az 

A 

Nr 


"B 


speed  of  light  (  3.00  x  108m/s) 

transmitted  wavelength  (m) 

effective  backscatter  cross  section  (m2) 

density  of  scatterers  (nr3) 

receiver  range  bin  length  (m) 

receiver  telescope  area  (m2) 

background  light  emission  rate  (photon/sec) 
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How  Much  Rayleigh  Signal  ? 


Consider  a  Rayleigh  lidar  system  transmitting  1  J  pulses  of  radiation  into 
the  sky  at  532  nm  and  receiving  the  backscattered  light  back  with  a  1-m 
telescope  at  10,  30,  60,  90  and  120  km. 

♦eff  10-32  m2 

p(z)  1025  m-3, 1024  m-3, 1022  m-3, 1020  m3, 1018  rrv3 


Az 

Efficiency 

Scattering  probability 
Receiver  probability 

Transmitted  Signal 

Photons  per  pulse 


1  km 
5% 

lO-3, 10-5, 10-7, 10-9, 10-11 

8x10-9,  9xl0-10,  2xl0-10, 10-10,  5xl0-n 


3xl018 


Expected  Signal 

~  12xl05  photons/pulse  at  10  km 
~  14xl02  photons/pulse  at  30  km  (1:860) 
~  3.0  photons/pulse  at  60  km  (1:470) 

~  15x10*3  photon/pulse  at  90  km  (1:200) 

~  75xl0*6  photon/pulse  at  120  km  (1:200) 


Temperature  profiles  measured  with  Rayleigh 
lidar  at  Poker  Flat  Research  Range. 
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How  Much  Resonance  Signal  ? 

Consider  a  Sodium  Resonance  lidar  system  transmitting  10  mJ  pulses 
of  radiation  into  the  sky  at  589  nm  and  receiving  the  backscattered 
light  back  with  a  1-m  telescope  at  90  km. 


p(z) 

Az 

Efficiency 

Scattering  probability 
Receiver  probability 


10-16  m2  (10-32  m2) 
1010  rrr3  (1020  m-3) 
1  km 
5% 

10-3  (10  9) 

10-io  (io-io) 


Transmitted  Signal 

Photons  per  pulse  3xl016  (3xl018) 


Expected  Signal 

~  150  (15xl0*3)  photon  received  per  pulse  at  90  k 


Sodium  Density  and  Potential  Temperature 
;d  by  Resonance  Lidar  at  Colorado 
rafeUniversity,  Fort  Collins,  CO,  USA. 


m 
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Why  Do  Signal  Levels  Matter  ? 

The  fundamental  uncertainty  in  the  measurement 

AN(z)  =  (N(z))1/2 

□arises  due  to  the  statistical  nature  of  the  detection  process  -  Poisson  Statistics. 

A  signal  level  of  100  photon  counts  has  an  inherent  uncertainty  of  10  photon 

counts  or  10%. 

□However,  if  a  signal  of  100  photon  counts  is  accompanied  by  a  background  of 
9900  photon  counts  yielding  a  total  uncertainty  of  (100+9900)1/2  =  100  and 

therefore  a  relative  uncertainty  of  100%. 

□To  increase  a  signal  quality  by  10  requires  we  increase  a  signal  magnitude  by  100. 

□We  can  play  with  Az  and  At  to  get  a  given  measurement  quality 

□  **  Resolution  vs.  Accuracy  ** 
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The  Ionosphere  and  Thermosphere 


MSISE-90 


100 


1000 


T emperature  (K) 


MSISE-90 


Density  (kg/m3) 


MSISE-90  MSISE-90  MSISE-90 


Altitude  (km) 


The  Aurorally-Modif  ied  Ionosphere 
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The  aurora  modifies  the 
composition  of  the  ionosphere. 
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Photo  courtesy  of  GI-UAF  by  Jan  Curtis 


W  a ve  n  u  m  be  r  (c  m 


The  Spectroscopy  of  Molecular  Nitrogen 


Molecular  Nitrogen  Term  Diagram 


> 


---i 

m 

<U 


PJ 


Partial  term  diagram  for 
molecular  Nitrogen. 
Transitions  of  interest 
include 

N2(A3S+U)  ->  N2(B3ng) 

the  "first  positive"  system 
N2(B3ng)  ->  N2(C3nu) 
the  "second  positive"  system, 

and 

N2+  (X2 1+g)  ->  N2+(B3  E+u) 
the  "first  negative"  system. 
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Choosing  an  Viable  Scatterer 

The  "first  negative"  transitions  (i.e.,  N2+  (X2I+g)  -»  N2+(B3E+U)) 
have  a  signif icantly  higher  Einstein  A-Coeff  icient  than  the  "first 
positive"  transitions  (i.e.,  N2(A3X+U)  N2(B3ng) ) 
and  the  "second  positive"  transitions  (i.e.,  N2(B3ng)  -»  N2(C3nu)). 

For  N2+,  ♦eff~5xl018  m2  (10  16  m2  for  Na) 

The  chosen  vibrational  transition  (■  =  0  ♦■  =  0  #■  =  0)  is  2.4 
times  more  optically  eff icient  than  the  next  most  efficient 
transition  (■  =  2  *  ■  =  3  *  ■  =  2). 

The  excited  oxygen  transition  0(3P)  -»  0(35°)  at  844.65  nm  has 
been  observed  but  the  lower  state  (i.e.,  0(3S°) )  is  very  short¬ 
lived  (~5  ns)  and  not  a  suitable  lower  state  for  lidar  observations  . 
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The  Spectroscopy  of  Molecules  vs.  Atoms 


Molecular  spectroscopy  has  vibrational  and  rotational  states. 
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Rotational  Population 


10' 


T  =  100  K 
J*  =  4 


\ . T  =  200  K 

J*  =  6 


l I  I  I  I  I  I  I  ll 


T  =  500  K 
J*  =  9 

i  i  i  i  i  L 


10  20  30  40 

Rotational  Quantum  Number,  J 


The  optimum  lidar 
wavelengths  vary  with 
rotational  population. 


Rotational  Population 


Temperature 

J 

Relative 

Population 

Transmit 

Wavelength 

Receive 

Wavelength 

100  K 

4 

14% 

390.781  nm 

391.428  nm 

200  K 

6 

10% 

390.604  nm 

391.486  nm 

500  K 

9 

6.4% 

390.303  nm 

391.537  nm 

1000  K 

13 

4.5% 

389.839  nm 

391.541  nm 
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Intensity  (Ry)  Intensity  (Ry)  Intensity  (Ry) 


Simulation  of  a  Case  Study 

6-7  February  1994 


Time  (LST  =  UT-9  h) 


A  meridian  scanning 
photometer  (MSP)  recorded 
the  intensity  of  the  skylight 
in  the  magnetic  meridian  on 
the  night  of  6-7  February 
1994  during  moderate  auroral 
activity. 

The  MSP  observations  allow 
determination  of  the  total 
energy  flux  and  characteristic 
energy  of  the  precipitating 
particles 
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Estimating  Concentrations  of  N2+ 


The  results  of  the 
MSP  analyses  are 
input  into  a  model  — 
that  simulates  the 
aurorally-modified  ^ 
atmosphere.  The 
model  estimates  the  H 
distribution  of  ^ 
molecular  nitrogen 
ions 


6-7  February  1994 


100 


300 

Half 


20  22  00  02  04  06 


Time  (LST) 
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A  Prototype  Lidar  System 


Proposed  N2+  XB  Lidar  System 


Transmitter 

Laser  Excimer-pumped  Dye 

Pulse  Energy  50  mJ 

Pulse  Repetition  Rate  20  pps 

Wavelength  390.303  nm 

Measurement  Resolution 
Temporal  Resolution  300  s 

Spatial  Resolution  5  km 


FAIRBANKS 


Receiver 

Telescope  Diameter  1  m 

Detector  Photomultiplier  Tube 

Filter  Bandwidth  0.3  nm 

Wavelength  391.537  nm 
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Lidar  Performance 


SNR 

50  40  30  20  10  0 


SNR 


Photon  Count 


The  expected  signal  varies  with 
auroral  activity.  The  best 
measurements  are  made  near  0057 
LST.  However,  the  lidar  system  is 
capable  of  making  measurements 
throughout  the  night. 


N  ‘  XB  Lidar  SNR 
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Altitude  (km) 


The  Advanced  Modular  Incoherent 
Scatter  Radar  (AMISR) 


Lidar  measurements  of  N2+  will  be 
nested  within  the  AMISR 
measurements  of  the  total  plasma. 
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Combined  Measurements 
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21SGS2.2-  21 3ftS3.S  UT 


25C 


'V 


100 


■i  ' 


300 


2S0  i 


214352.0  -  214*54.2  UT 


ISO 


south 


N0PTH 


21 4723.2  -  21  5024,4  UT 


50 

oL 


21 4022.5-  214323.9  UT 


-100  0 


tOO  -100  0  100 

GROUND  DISTANCE  |Km) 


For  ionization  densities  of  4x105  cm-3  at  115  km  AMISR  will  make 
10%  measurements  at  3  km  resolution  every  2  seconds  while  the 
N2+  lidar  will  make  10%  measurements  at  1  km  resolution  every  30 

seconds. 
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INTEGRATION  TIME  (s) 


E-Region  Chemical  Assimilation 
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Beyond  Concentration  Profiling 


Rotational  State  Resonance  Lidar 


LIDAR  TRANSMITTER 


LIDAR  RECEIVER 


A  dual  laser  lidar  system  employing  solid-state  lasers  could  be  used  to  prof  ile  two 
rotational  states  simultaneously  and  hence  study  the  energy  deposition  in  the 

auroral  ionosphere. 
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Conclusions  and  Discussion 


N2+  lidar  appears  feasible  for  studies  of  the  auroral  thermosphere. 

Range-resolved  measurement  of  a  specif  ic  ion  extends  the  scope  of 
current  radar/photometer  experiments. 

The  lidar  differs  in  several  technical  aspects  from  "traditional" 
atomic  resomnce  lidars. 

Operational  techniques  will  be  developed  to  yield  robust  lidar 
measurements. 

Assimilation  techniques  will  be  developed  to  yield  robust  retrievals 
of  ionospheric  species. 
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